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ABSTRACT 


Graphical  computer  simulations  provide  a  means  through  which  weapon  prototyping 
and  tactical  evaluations  can  be  conducted  at  low  cost,  without  the  risks  associated  with  the 
movement  of  equipment  and  firing  of  weapons.  Because  of  the  widespread  use  of  ballistic 
munitions  in  the  armed  forces,  a  fundamental  aspect  of  the  implementation  of  such  mUitary 
simulations  is  a  physical  model  that  governs  ballistic  behavior.  The  modified  point-mass 
trajectory  model  is  used  to  implement  ballistic  trajectories  within  NPSNET,  a  real-time, 
graphical,  three-dimensional  simulation.  A  parallel  algorithm  is  used  to  simulate  the  visual 
characteristics  of  shrapnel-producing  explosions. 

A  special  case  of  ballistic  trajectories  involves  the  application  of  indirect  fires.  When 
a  projectile  travels  along  a  curved  path  to  the  target  area,  rather  than  being  propelled 
directly  along  the  line  of  sight,  much  greater  ranges  can  be  achieved.  This  makes  it  possible 
to  fire  upon  an  enemy  without  directly  exposing  the  firing  elements  to  harm.  As  a  result  of 
these  increased  ranges,  it  is  generally  not  possible  for  the  firing  element  to  acquire  its  own 
targets.  Thus,  an  additional  player  is  required  to  represent  this  tactic  in  a  virtual  world:  the 
forward  observer.  An  expert  system  is  presented  that  mimics  the  cognitive  contributions 
of  a  human  forward  observer. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Public  pressure  to  reduce  expenditures  for  national  defense,  coupled  with  a  desire  to 
field  state  of  the  art  defense  technology  in  a  timely  manner  have  led  the  Department  of 
Defense  to  examine  ways  to  “short-circuit”  the  traditionally  time-consuming  process  of 
designing  and  fielding  new  weapon  systems  (Atwood,  1991).  The  ideal  technique  would 
allow  for  rapid  prototyping  and  testing,  and  permit  changes  to  original  specifications 
without  need  of  an  extensive  re-engineering  effon. 

B.  SYNTHETIC  ENVIRONMENTS 

1.  General. 

A  rapidly  maturing  technology  that  could  potentially  produce  significant  savings 
in  the  pursuit  of  low-cost  weapons  prototyping  is  the  development  of  graphical  computer 
simulations.  Variously  known  in  the  literature  as  “artificial  reality”,  “virtual  worlds”, 
“synthetic  environments”,  and  “virtual  reality”,  the  science  of  simulation  has  advanced 
dramatically  since  the  first  purely  statistical  combat  models.  The  ability  to  create  an 
electronically  manipulable  milieu  inside  which  visible  representations  of  objects  can  be 
made  to  imitate  the  real-world  objects  that  they  resemble  has  achieved  some  notoriety  in 
the  lay  press  of  late,  but  to  be  sure,  some  of  what  sounded  like  shades  of  Flash  Gordon  as 
few  as  five  years  ago  has  come  to  pass. 

2.  Flight  simulators. 

Graphical  simulations  have  already  seen  extensive  application  in  some  areas  of 
the  military  genre.  The  earliest  examples  of  these  included  aircraft  flight  simulators  used  to 
train  pilots  prior  to,  or  even  in  lieu  of,  actual  aircraft  use.  A  mock-up  of  an  aircraft  cockpit 
was  slaved  to  a  computer  that  presented  the  operator  with  a  display  that  changed  the 
viewpoint  and  orientation  according  to  inputs  received  from  the  controls.  Updating  the 
display  at  a  sufficiently  rapid  rate  produced  the  illusion  of  motion  through  a  three- 
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dimensional  space  according  to  the  rules  of  flight.  This  first  use  of  graphical  simulations 
remains  a  prevailing  use  today,  albeit  in  a  much  more  sophisticated  form.  Zyda  et  al  have 
developed  an  inexpensive  simulator  for  the  purpose  of  modeling  the  flight  of  a  guided 
munition,  and  testing  the  application  of  this  weapon  to  different  tactical  scenarios  (Zyda, 
McGhee,  Smith,  and  Streyle,  1987,  pp.  10-11). 

3.  Two-dimensional  ground  combat  models. 

The  first  uses  of  graphical  simulations  for  the  purpose  of  modeling  the  actions  of 
ground  forces  (many  of  which  remain  in  use)  used  the  computer  display  as  a  two- 
dimensional  automated  wargaming  platform.  The  computer  accepted  input  in  textual  form 
or  from  a  mouse  that  corresponded  to  actions  in  the  simulation.  The  display  was  used 
mainly  to  keep  track  of  units’  positions  from  an  overhead  (or  “god’s  eye”)  view  of  the 
simulated  terrain.  Engagements  between  ground  units  were  modeled  by  some  stochastic 
means  based  on  the  empirical  characteristics  of  the  weapons  involved.  Examples  of  these 
kinds  of  simulators  include  ARTBASS  (Department  of  the  Army,  1987),  and  JANUS(A) 
(Department  of  the  Aamy,  1986). 

A  disadvantage  of  these  two-dimensional  programs  is  that  they  offer  a  very  large- 
scale  view  of  the  virtual  world  at  hand.  The  real  world  is  of  course  decidedly  three- 
dimensional.  A  single  user  of  such  a  system  therefore  is  typically  not  presented  with 
sufficient  stimuli  to  persuade  him  that  he  is  to  any  extent  participating  in  the  simulation  that 
the  computer  is  modeling.  On  the  other  hand,  when  the  objective  of  the  sunulation  is  to  train 
high-level  staffs,  the  two-dimensional  representation  is  a  very  familiar  one,  and  thus 
completely  satisfactory. 

4.  Three-dimensional  models. 

State  of  the  art  military  graphical  simulations  are  real-time,  three-dimensional, 
interactive  systems.  They  include  the  capability  to  visualize  world  objects  such  as  teiram, 
vehicles,  buildings,  and  vegetation  in  depth.  The  user  can  move  over  the  terrain  and  look 
around.  As  the  user’s  position  and  orientation  changes,  the  display  device  changes 
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accordingly,  presenting  a  view  that  corresponds  approximately  to  the  scene  the  user  would 
obtain  when  similarly  oriented  in  an  authentic  environment  of  the  same  construction.  Since 
the  display  is  updated  in  real  time,  the  user  has  the  sensation  of  actually  experiencing  the 
events  that  are  being  simulated.  NPSNET  (Zyda  &  Pratt,  1992),  SIMNET  (Pope,  1989), 
and  JANUS-3D  (Walter  &  Warren,  1992)  are  all  simulators  of  this  category. 

The  three-dimensional  nature  of  the  latter  category  of  simulations  gives  them  a 
powerful  advantage  over  their  two-dimensional  counterparts.  The  3D  simulation,  being  the 
more  general,  has  the  greatest  efficiency  in  terms  of  utilization  of  the  optical  bandwidth. 
There  are  many  cues  within  a  3D  context  that  are  accessible  for  encoding  information  about 
the  simulated  environment  that  are  unavailable  or  ludicrous  in  2D.  For  example,  shadows 
can  offer  hints  to  an  observer  about  an  object’s  spatial  orientation  (Chamiak  &  McDermott, 
1985,  pp.  129-130).  The  increased  capacity  of  information  transmission  allows  for  the 
representation  of  more  complex  objects,  and  in  greater  physical  detail.  This  level  of 
specificity  is  essential  to  the  creation  of  believable  synthetic  environments. 

5.  Benefits. 

The  versatility  of  the  typical  3D  simulation  makes  it  especially  well-suited  for  a 
variety  of  cost-saving  military  applications.  Flight  simulations  are  a  classic  example.  In  a 
networked  system,  the  simultaneous  interaction  of  large  numbers  of  individuals  is  possible. 
The  availability  of  world-wide  telecommunications  networks  makes  it  viable  for  units  on 
one  side  of  the  globe  to  engage  in  mock  combat  against  units  on  the  other  side.  Through 
use  of  autonomous  forces,  troops  may  train  for  missions  involving  combat  with  an  enemy 
without  need  for  a  human  opposing  force’s  participation. 

Aside  from  the  obvious  desirability  of  enhanced  cost-effectiveness,  a  far  more 
attractive  characteristic  emerges  from  the  consideration  of  graphical  simulations:  improved 
safety.  The  act  of  training  human  beings  to  conduct  armed  warfare  is  inherently  very 
dangerous.  It  invariably  requires  frequent  interaction  at  close  quarters  with  explosives, 
heavy  equipment,  and  petroleum  products.  Training  of  any  significant  duration  produces 
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fatigue,  which  can  lead  to  decreased  motor  capacity  and  poor  judgement.  Exposure  to  the 
weather  may  further  degrade  a  fighting  man’s  ability  to  take  actions  to  protect  himself 
against  unforeseen  hazards.  If  training  can  be  conducted  in  a  virtual  setting,  these  threats  to 
life  and  limb  can  be  eliminated.  Individuals  can  be  trained  in  the  operation  or  maintenance 
of  various  types  of  equipment  without  risk  to  themselves  or  the  machinery.  Mortal  combat 
can  be  simulated  without  actual  bloodshed,  either  accidental  or  intended.  The  3D 
interactive  virtual  world  offers  the  greatest  hope  among  diverse  technologies  for 
maintaining  a  credible  defensive  capability  in  the  face  of  waning  public  support  for 
adequate  funding. 

C.  AGENTS  OF  FORCE  IN  A  VIRTUAL  WORLD 

In  most  simulations,  it  is  useful  to  monitor  the  participant’s  status,  and  render  a  report 
at  certain  intervals  with  some  kind  of  comparative  analysis.  In  game-type  simulations,  this 
is  known  as  “keeping  score”.  For  the  military  force-on-force  analogue,  a  means  is  needed 
to  determine  the  unit  or  individual  that  will  prevail  in  the  course  of  an  engagement  In  real 
combat,  of  course,  the  force  that  is  able  to  inflict  the  greatest  amount  of  damage  upon  its 
opponent  is  the  victor.  Thus,  any  simulation  that  models  the  interaction  of  military  units 
necessarily  implies  the  implementation  of  various  kinds  of  weapons,  since  they  are  the 
principal  means  by  which  deadly  force  is  applied. 

The  agents  of  force  presently  in  use  can  be  categorized  according  to  the  characteristics 
of  their  flight  through  the  air.  Ballistic  munitions  travel  along  a  roughly  parabolic 
trajectory,  with  only  naturally  occurring  forces  such  as  gravity,  air  friction,  and  wind 
determining  their  course  after  launch.  Guided  and  rocket-propelled  munitions  rely  upon  the 
action  of  man-made  forces  such  as  the  movement  of  flight  control  surfaces,  or  the  forces 
generated  by  the  expulsion  of  gases  from  a  propellant  source  during  the  course  of  their 
flight  to  determine  heir  point  of  impact.  Both  kinds  of  devices  have  found  widespread  use 
in  the  armed  forces  of  the  world.  Weapons  of  the  ballistic  variety,  however,  are  by  far  the 
most  prevalent  among  ground  forces.  Main  battle  tanks,  infantry  fighting  vehicles,  most 
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field  artillery,  mortars,  crew-served  machineguns,  and  of  course,  the  infantryman’s  rifle  are 
all  examples  of  ballistic  weapons. 

The  focus  of  this  thesis  is  an  implementation  of  a  ballistic  weapons  model  within  the 
context  of  a  synthetic  environment.  Techniques  appropriate  for  the  simulation  of  in-flight 
and  terminal  characteristics  are  developed,  as  well  as  support  modules  necessary  for  the 
special  case  of  indirect  ballistic  systems.  The  most  significant  characteristic  of  these 
techniques  are  that  they  operate  in  rci-l-time,  providing  feedback  to  the  event  control  loop 
with  sufficient  rapidity  to  allow  the  perception  of  cause  and  effect  necessary  to  maintain 
the  persuasiveness  of  the  simulation. 
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II.  NPSNET  AND  OTHER  PREVIOUS  WORK 

A.  NPSNET 

NPSNET  is  a  real-time,  graphical  simulation.  Its  focus  is  the  interaction  of  ground  and 
air  vehicles  in  a  combat  environment  (Zyda  &  Pratt,  1992).  Any  number  of  humans  can 
participate  through  use  of  networked  graphics  workstations.  At  present,  NPSNET  provides 
for  a  company  of  autonomous  vehicles  that  can  be  assigned  missions  in  the  manner  of 
aggressor  forces  (Culpepper,  1992).  By  default,  the  remainder  of  the  vehicles  that  exist  in 
this  virtual  world  behave  in  an  autonomous  fashion,  albeit  statically.  These  vehicles  are 
assigned  random  velocities  and  orientations  upon  program  start-up,  and  move  over  the 
simulated  terrain  accordingly.  If  fired  upon,  they  will  alternatively  return  fire,  or  retire  at 
high  speed.  The  capability  exists  to  record  and  playback  sequences  of  interaction,  and  to 
view  these  sequences  from  any  viewpoint  available  in  the  world  space. 

The  primary  means  by  which  deadly  force  is  applied  in  NPSNET  is  through  direct, 
line-of-sight  fires.  When  a  player  sees  an  enemy  vehicle  on  the  screen,  he  visually  orients 
himself  so  as  to  be  able  to  fire  in  the  direction  of  the  target.  A  projectile  object  is  generated 
at  the  moment  of  firing  which  follows  a  trajectory  that  is  linear  by  default.  Collision 
detection  of  the  projectile  is  performed  along  its  path,  and  a  vehicle  is  “killed”  if  the 
projectile  passes  within  some  specified  distance. 

NPSNET  provides  an  excellent  platform  upon  which  to  implement  realistic  models  of 
ballistic  weapons.  The  availability  of  visual  representations  of  a  plethora  of  vehicles  and 
objects  sets  the  stage  for  a  virtual  firing  range.  Numerical  models  of  the  flight  of  various 
kinds  of  projectiles  can  be  constructed,  and  their  efficacy  tested  by  direct  observation 
within  the  simulator.  The  behavior  of  entities  intended  to  control  the  application  of  direct 
and  indirect  fires  may  be  examined  visually  as  well. 

B.  SURVEY  OF  PREVIOUS  WORK 

The  computation  of  ballistic  trajectories  using  digital  computers  has  a  long  and 
distinguished  history.  The  venerable  ENIAC  was  designed  specifically  for  the  purpose  of 
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automating  the  tedious  process  of  calculating  weapons  data  for  the  United  States  Army 
(Tanenbaum,  1990,  p.  16).  Attempts  to  model  ballistic  weapons  within  the  context  of  a  real¬ 
time  graphical  simulation,  however,  have  begun  fairly  recently. 

1.  Graphical  Simulators. 

a,  SIMNET 

SIMNET  (Garvey  &  Monday,  1988)  is  a  sophisticated  three-dimensional 
military  simulator.  It  makes  use  of  specially  designed  hardware  and  software  to  provide  a 
networked  interactive  virtual  world.  The  flight  of  ballistic  munitions  are  modeled  using  a 
modified  static  linear  algorithm.  Terminal  visual  effects  are  obtained  by  the  intermittent 
presentation  of  texture  maps.  Kills  are  determined  by  a  stochastic  selection,  using  historical 
distributions  of  targeting  probable  errors.  Direct  fires  are  controlled  by  either  of  the 
following  methods: 

•  by  console  operators  playing  the  role  of  a  vehicle  commander, 

•  by  Semi-Autonomous  Forces  (SAF)  entities, 

•  by  human  SAF  commanders. 

Indirect  fires  in  SIMNET  must  be  controlled  by  the  SAF  commander,  a 
human  console  operator  that  supervises  the  interaction  of  SAF  forces  with  forces  controlled 
by  other  humans  using  the  system. 

b.  NPSNET 

In  previous  versions  of  NPSNET,  the  paths  of  ballistic  projectiles  were 
calculated  by  direct  interpolation  between  the  launch  position  and  the  point  of  impact. 
Thus,  indirect  fires  were  not  modeled.  NPSNET  also  makes  use  of  an  animated  sequence 
of  texture  maps  to  give  a  visual  representation  of  terminal  weapons  effects.  Kills  are 
determined  by  performing  3D  collision  detection  within  the  boundaries  of  two  invisible 
concentric  polygons  surrounding  the  projectile.  Objects  detected  as  being  within  the 
outermost  polygon  are  assessed  damage;  objects  within  the  inner  polygon  are  destroyed 
(Osborne,  1991).  Monahan  has  developed  a  technique  by  which  the  kinematics  of  particles 
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generated  in  the  decomposition  of  an  exploding  target  can  be  modeled  based  on  the  force 
environment  in  which  the  explosion  takes  place  (Monahan,  1991).  This  technique  proved 
to  be  too  computationally  intensive  to  implement  within  NPSNET,  given  the  state  of  the  art 
in  graphics  hardware.  Branley  and  Culpepper  devised  the  means  to  create  instances  of 
platoon-sized  elements  in  NPSNET  which  could  react  autonomously  to  the  presence  of 
other  entities  in  the  virtual  world  according  to  previously  established  mission  instructions 
(Branley,  1992)  (Culpepper,  1992). 

c.  POST 

The  Forward  Observer  Simulation  Trainer  (POST)  was  created  by 
Drummond  and  Nizolak  (Drummond  and  Nizolak,  1989).  It  is  a  graphical  application  that 
simulates  the  interface  of  the  Digital  Message  Device,  which  is  the  means  by  which  Army 
indirect  fire  spotters  communicate  with  the  firing  elements.  Since  the  primary  intent  of  their 
work  was  to  develop  an  automated  means  for  soldiers  to  train  with  this  particular  interface 
device,  the  representation  of  the  projectile’s  actual  trajectory  was  not  developed;  rather,  the 
projectiles  appeared  at  the  location  specified  by  the  observer  at  the  appropriate  time. 

2.  Particle  systems  of  ballistic  entities. 

Reeves  used  a  particle-based  approach  to  model  the  terminal  visual  effects  of  an 
exploding  “wall  of  fire”  (Reeves,  1983,  p.  365).  This  was  a  successful  implementation 
using  the  abstraction  of  a  dynamic  particle  system  for  a  graphical  simulation.  Later,  he 
extended  the  technique  to  include  objects  that  reflect  light,  as  well  as  emit  it  (Reeves  and 
Blau,  1985  p.  313).  More  recently,  Loke  et  al  used  a  similar  technique  for  representing  the 
visual  characteristics  of  fireworks  (Loke,  Tan,  Seah,  and  Er,  1992,  p.  33). 
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III.  EXTERIOR  BALLISTICS 


A.  DEFINITION 

The  study  of  exterior  ballistics  pertains  to  the  description  of  the  behavior  of  a  ballistic 
projectile  as  it  flies  through  the  air.  More  specifically,  it  is  concerned  with  the  period  of 
time  between  the  moment  the  projectile  leaves  the  barrel  of  the  weapon  from  which  it  was 
fired,  and  the  moment  that  it  explodes  in  the  air  or  strikes  the  ground.  In  a  virtual  world, 
being  able  to  establish  the  position  of  a  projectile  as  a  function  of  time  is  fundamental  to 
the  successful  modeling  of  any  sort  of  weapon  system.  This  capability  is  necessary  in  order 
to  render  a  scene  that  contains  the  projectile,  as  well  as  to  perfonn  collision  detection  with 
the  other  objects  in  the  world.  In  order  to  fully  understand  what  happens  following  the 
departure  of  the  projectile  from  the  weapon,  a  brief  introduction  to  interior  ballistics  is  in 
order.  Figure  1  shows  the  general  characteristics  of  an  artUlery  projectile. 


Figure  1;  Separate-Loading  Artillery  Projectile  (Shown  Without  Fuze). 

Interior  ballistics  is  the  study  of  the  processes  that  take  place  uiside  the  bore  of  a 
weapon.  In  general,  these  processes  are  responsible  for  placuig  the  projectile  in  motion,  and 
to  a  large  extent  determine  the  overall  behavior  of  the  projectile  during  flight.  The  majority 
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of  ground-based  weapons  in  the  United  States  anned  forces  are  spin-stabilized.  That  is, 
their  projectiles  revolve  rapidly  about  their  longitudinal  axis  of  symmetry  while  in  flight. 
Tlie  spinning  of  the  projectile  provides  stability  to  the  trajectory,  as  a  result  of  tlie 
gyroscopic  effect.  The  spin  is  developed  in  the  following  way:  when  the  projectile  is 
rammed  into  the  breech,  the  rotating  band  (in  the  case  of  most  artillery  shells)  or  the  widest 
part  of  the  round  (in  the  case  of  munitions  lackmg  a  rotating  band)  contacts  the  lands  and 
grooves  of  the  bore.  When  the  weapon  is  fired,  the  force  developed  by  the  rapidly 
expanding  gases  of  the  burning  propellant  push  the  projectile  forward,  causing  the  rotating 
band  to  be  engraved  by  the  rifling  of  the  tube.  As  the  projectile  travels  the  length  of  the  tube, 
spin  is  imparted  as  it  twists  over  the  lands  and  grooves.  Thus  when  the  round  exits  the  tube, 
its  motion  has  both  translation  and  rotation  characteristics. 

B.  NEWTONIAN  MECHANICS 

The  discipline  of  classical  mechanics  provides  us  with  the  simplest  method  of 
modeling  the  motion  of  a  projectile  through  space.  The  term  “projectile  motion”  is  usually 
applied  to  the  motion  of  an  object  acting  solely  under  the  influence  of  the  force  of  gravity. 
The  position  of  an  object  that  moves  in  this  fashion  may  be  described  as  a  function  of  time 
in  a  Cartesian  coordinate  system  by  the  following  well-known  formulae; 


(Eq3.1) 

v(0  =  +  + 

(Eq3.2) 

1  2 

^(0  = 

^  *-0  z  *■ 

(Eq3.3) 

Since  we  assume  that  the  only  force  acting  upon  the  projectile  is  the  force  due  to 
gravity,  we  can  then  conclude  that  =  a,  =  0,  and  that  a^.  =  g.  where  g  is  the  standard 
gravitational  acceleration  of  the  Earth.  Equations  3.1.  3.2,  and  3.3  therefore  reduce  to; 
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(Eq  3.4) 


Xit)  + 

y(t)  =  }'o  +  (Eq  3.5) 

z(0  =  Zq  +  v^^/  (Eq3.6) 

The  fact  that  these  equations  have  only  one  independent  variable,  namely  t,  implies 
that  the  position  function  can  be  implemented  in  a  computer  program  with  a  minimum  of 
parameter  passing.  This  means  that  such  a  function  can  be  called  many  times  throughout 
the  course  of  a  computer  graphics  program’s  screen  refresh  cycle  without  greatly  increasing 
the  program  overhead.  For  these  reasons,  it  is  clear  that  Equations  3.4  -  3.6  are  very  well- 
suited  for  use  in  a  real-time  simulation  requiring  minimal  delays  in  calculating  object 
position. 

While  these  equations  of  projectile  motion  have  the  desirable  characteristics  of  being 
straightforward  and  rapidly  computable,  they  are  highly  simplified  representations  of  a 
complex  phenomenon.  As  such,  their  use  is  limited  in  applications  which  require  a  high 
degree  of  correspondence  between  physical  reality  and  the  simulated  environment,  such  as 
weapons  prototyping  or  tactics  evaluation.  Table  2  shows  a  comparison  of  values  obtained 
using  this  technique  and  values  from  U.S.  Army  artillery  firing  tables  (Department  of  the 
Army,  1983). 


TABLE  1:  SIMPLIFIED  MODEL  DATA  VS.  FIRING  TABLE  DATA 


Charge^ 

Quadrant 

elevation 

(mils) 

Observed 

impact 

(deflection, 

range) 

Firing  table 
(deflection, 
range) 

CEP 

'y 

(meters) 

Error 

(number  of 
CEPs) 

3W 

266.67 

(4,5513) 

(22, 4056) 

10.66 

136.7 

3W 

533.33 

(6,  8794) 

(84.5,  6588) 

17.49 

126.2 

3W 

800.0 

(7,  10000) 

(207,  7300) 

21.35 

126.8 
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TABLE  1;  SIMPLIHED  MODEL  DATA  VS.  FIRENG  TABLE  DATA  (CONTINUED) 


Charge^ 

Quadrant 

elevation 

(mils) 

Observed 

impact 

(deflection, 

range) 

Firing  table 
(deflection, 
range ) 

CEP 

(meters)' 

Error 

(number  of 
CEPs) 

7W 

266.67 

(5,  17490) 

(69,  8983) 

9.8 

868.3 

7W 

533.33 

(9, 29522) 

(220, 

13154) 

14.97 

1093.5 

7W 

800.0 

(10,33127) 

(441, 

14700) 

16.73 

1101.7 

1.  PropeUant  charge  is  identified  by  a  zone  (a  numeral  from  one  to  eight,  indicating  the  general 
magnitude  of  the  charge),  and  a  mode!  (a  particular  size  and  shape  of  propellant  grain).  In  this 
case,  the  model  is  M4A2  white  bag. 

2.  Circular  Error  Probable. 

C.  MODIFIED  POINT-MASS  MODEL 

The  assumption  that  no  forces  other  than  gravity  act  upon  a  projectile  in  free  flight  was 
empirically  shown  to  be  an  overly  broad  simplification.  The  fact  is  that  several  other  forces 
do  act  on  spin-stabilized  projectiles  tliroughout  their  trajectory,  namely: 

1.  Air  drag. 

Projectiles  are  not  perfectly  smooth.  Defects  in  the  casting  of  the  projectile 
casing,  gouges  in  the  area  of  the  rotating  band  caused  by  contact  with  the  forcing  cone  of 
the  bore,  and  variations  in  the  formation  of  the  molecular  matrix  of  the  projectile  skin  are 
all  perturbations  which  impede  the  flow  of  air.  As  a  result,  the  motion  of  the  projectile 
through  the  air  generates  a  shearing  force  as  the  surface  of  the  round  slides  past  the  body 
of  air  that  surrounds  it. 

2.  Lift. 

The  angular  deflection  in  the  vertical  plane  between  the  longitudinal  axis  of  the 
projectile  as  it  leaves  the  barrel  and  an  imaginary  line  that  passes  through  the  base  of  the 
trajectory  is  called  the  quadrant  elevation.  At  point  of  launch,  the  forward  portion  of  the 
round's  surface  area  is  presented  such  that  as  it  moves  through  the  fluid  medium,  the  air 
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moves  at  differing  velocities  over  the  upper  and  lower  portions  of  the  projectile.  Tliis 
difference  in  air  velocity  creates  a  pressure  differential  that  imparts  a  corresponding  force 
to  the  projectile,  known  as  dynamic  lift.  Increasing  quadrant  elevation  results  in  an  increase 
in  lift,  in  much  the  same  way  as  increasing  the  angle  of  attack  of  an  airfoU  does. 

3.  Equilibrium  yaw. 

As  the  projectile  negotiates  the  rifling  of  the  bore  in  its  travel  during  launch,  the 
bourrelet  of  the  round  rides  over  the  lands  and  grooves.  Nicks  and  pits  in  the  lands,  as  well 
as  slight  imperfections  in  the  machining  of  the  bourrelet  exert  unpulsive  forces  whicli  cause 
the  round  to  precess  as  it  spins.  This  precession  eventually  damps  out.  but  results  in  a 
permanent  deflection  from  the  projectile’s  original  orientation.  This  difference  produces  a 
lateral  force  in  a  fashion  similar  to  the  way  that  lift  exerts  a  venical  force  on  the  round.  The 
change  in  the  trajectory  resulting  from  this  lateral  force  is  called  drift. 

4.  Coriolis  effect. 

Many  conditions  used  in  the  development  of  equations  of  motion  for  freely 
failing  bodies  assume  fairly  short  times  of  flight.  A  classic  example  is  that  of  a  baseball. 
However,  in  the  case  of  military  ammunition,  ra.nges  are  typically  very  large,  sometunes  on 
the  order  of  tens  of  kilometers.  Depending  on  the  quadrant  elevation,  times  of  flight  can  be 
as  long  as  a  minute  or  more.  As  a  result,  the  rotation  of  the  earth  has  an  effect  on  tlie 
perceived  trajectory,  in  that  the  time  differential  between  launch  and  landing  is  large 
enough  to  make  the  radial  motion  of  the  earth  noticeable.  The  effect  of  this  rotation  is  to 
cause  an  increase  or  decrease  in  the  observed  achieved  range,  depending  on  the  latitude  of 
the  weapon. 

The  point  mass  trajectory  model  was  developed  by  Lieske  and  Reiter  to  provide  a 
means  of  calculating  artillery  trajectories  having  a  computational  burden  less  than  the 
complete  rigid  body  model,  yet  stiU  incorporating  the  effects  of  the  forces  described  above. 
(Lieske  &  Reiter,  1966,  p.  19).  The  modified  point-mass  model  was  later  developed  to 
incorporate  the  effects  of  thrust  from  rocket-assisted  projectiles. 
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According  to  the  modified  point  mass  model  (Lieske.  1973,  pp.  7-8),  the  equation  of 
motion  for  a  fi^ee-flight  projectile  is  given  by; 


du. 

dt 


rpHm^v._i 


n~ 

%Cm 


+ 


^pKd^Lv^i_  ^ 
^  8m 


+  (x^ja^  +  g  +  A-i- 


^-pKd^pQC^^  ^ 

^  o. 

8  m 


(Eq  3.7) 


V  J 

where  U  is  the  velocity  of  the  projectile  with  respect  to  the  ground,  v  is  the  velocity  of  the 
projectile  with  respect  to  the  air.  and  Ot^  is  the  estimate  of  the  yaw  of  repose  at  the  ( i-1  )th 
time  step.  See  Appendix  A  for  a  complete  listing  of  the  meanings  of  the  remaining 
variables. 

Since  this  algorithm  is  intended  for  use  in  a  real-time,  graphics-intensive  application, 
the  following  assumptions  are  made  to  simplify  Equation  3.7; 

•  Mass  of  the  projectile  is  constant  ( m,  =  m). 

•  Coriolis  effect  is  negligible  ( A  =  (0,0,0)). 

•  The  force  of  gravity  acts  only  along  an  axis  that  is  perpendicular  to  the  trajectory. 

•  No  winds  are  present  (U  =  v ). 

These  assumptions  yield  the  following  equation  of  motion; 


dt 


prtv.-i 


8C 


+ 


'pnd^Lvj_  j 
8  m 
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-pndpQC^ 


(%xv,._i) 


(Eq3.8) 


Using  these  same  assumptions,  is  given  by: 


-8/p  d\- 


(Eq3.9) 


p,  which  is  the  magnitude  of  tlie  velocity  with  which  the  projectile  rotates  about  its 
longitudinal  axis  is  obtained  by  integrating; 


^  = 
dt 


pKct^  pc  I 


lEq  3.10) 


The  initial  spin  velocity,  pQ  is  given  by; 


Po= 


(Eq3.11) 


Equation  3.10  is  a  linear  first-order  differential  equation  in  t,  which  may  be  solved 
readily.  For  the  sake  of  notational  simplicity,  let 


(Eq3.12) 


The  general  solution  then  is  given  by; 

dr 

p{t)  -  ce  ^ 


p{t)  =  ce 


(Eq3.13) 


(Eq3.l4) 


Substituting  the  initial  condition  p(0)  =  V(j  gives  the  particular  solution,  namely: 
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P(0  = 


2KVQe 

~~bd 


at 


(Eq3,15) 


D.  IMPLEMENTATION 

An  object-oriented  approach  is  used  in  the  implementation  of  the  modified  point  mass 
model.  This  allows  for  an  intuitive  understanding  of  the  code  used  in  the  simulation.  A 
summary  of  the  major  object  classes  follows. 

1.  Class  WorldObject. 

This  is  the  base  class  for  objects  that  are  intended  to  have  some  physical 
representation  in  the  virtual  world.  As  a  result  of  the  pure  virtual  member  function  draw  { ), 
it  is  an  abstract  class,  and  thus  can  only  be  u.sed  to  derive  additional  subclasses.  Very’ 
rudimentary  levels  of  functionality  are  provided,  including  the  ability  to  set  and  report 
object  location  and  orientation,  as  well  as  a  tagging  scheme  to  uniquely  identify  objects 
without  need  of  using  an  enumerated  type. 

2.  Class  Particle. 

Originally  designed  for  use  in  modeluig  particle  systems,  this  class  turned  out  to 
be  equally  suitable  in  representing  projectile  objects.  It  is  derived  from  class  WorldObject, 
and  is  polymorphic  by  the  addition  of  data  elements  pertaining  to  the  object’s  physical 
characteristics,  such  as  initial  velocity  and  lifetime,  and  by  its  overriding  definition  of 
drawi).  The  definition  of  this  function  determines  the  visual  characteristics  of  the  particle. 
By  default,  objects  of  type  Particle  are  displayed  using  the  Naval  Postgraduate  School 
Object  File  Format  (NPSOFF)  utilities.  See  (Zyda.  Wilson,  Pratt,  and  Monahan,  1991 )  for 
a  detailed  development  of  the  capabilities  of  NPSOFF,  and  (Wilson,  September  1992)  for 
a  powerful  superset  of  object-oriented  extensions  to  the  original  design.  Tire  remaining 
attributes  of  this  class  are  discussed  in  detail  in  Chapter  IV. 
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3.  Class  Ballistic_co€fficient. 

The  equations  employed  by  the  modified  point  mass  model  make  use  of  a  number 
of  dimensionless  power  products  that  each  capture  the  state  of  some  physical  characteristic 
of  the  projectile’s  flight  at  a  given  time.  All  were  detennined  experimentally,  and  most  are 
functions  of  mach  number.  Some  are  also  dependent  upon  the  quadrant  elevation.  The 
intervals  of  mach  number  over  which  each  aerodynamic  attribute  is  defined  vary  for  each 
value,  and  an  updated  value  is  needed  for  each  iteration  of  the  main  equation  of  motion.  It 
is  therefore  convenient  to  group  sets  of  these  parameters  that  are  applicable  to  a  given  mach 
number  into  a  single  data  structure.  Objects  of  class  BalIistic_coefficient  contain  storage 
for  each  of  the  following  attributes; 

•  Yaw  drag. 

•  Drag  force. 

•  Lift  force. 

•  Overturning  moment. 

•  Magnus  force. 

•  Spin  damping  moment. 

•  Yaw  lift  force. 

•  Ballistic  coefficient  for  reference  mass. 

•  Lift  factor. 

The  last  two  items  are  not  coefficients  which  characterize  a  ballistic  trajectory; 
rather,  they  capture  various  physical  aspects  of  the  projectile  itself  that  relate  to  its  fonn  and 
mass  distribution. 

4.  Class  Physical_coefficlent. 

Because  the  intervals  over  which  the  coefficients  that  combine  to  form  a 
particular  aerodynamic  product  are  variable  in  size,  an  adaptive  technique  is  necessary  to 
extract  the  appropriate  value,  given  the  entry  parameter  of  mach  number.  Class 
Physical_coefficient  provides  a  means  of  encapsulating  one  interval  of  coefficients  and  the 
ceiling  mach  number  for  which  that  set  is  valid.  A  detailed  explanation  of  the  method  used 
to  calculate  the  products  will  be  given  in  Section  6. 
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5.  Data  structures. 


a.  Static  parameters 

The  ballistic  coefficient  and  the  lift  factor  used  in  the  equations  of  motion  lue 
functions  of  the  propellant  charge  and  quadrant  elevation,  and  are  independent  of  its 
velocity.  Since  these  parameters  do  not  vary  over  the  flight  of  the  projectile,  they  can  be 
calculated  only  once  at  the  beginning  of  the  integration,  and  stored.  The  entry  aigument 
used  to  determine  the  coefficients  needed  to  calculate  these  values  is  propellant  charge.  For 
separate-loading  artillery  ammunition,  charge  is  characterized  by  a  number  from  one  to 
eight,  and  a  color,  either  green  or  white.  Because  each  combination  of  these  attributes 
produces  a  discrete  value,  charge  is  a  more  convenient  value  to  use  as  a  table  index  than 
muzzle  velocity,  which  is  a  floating-point  value.  Three  lookup  tables  indexed  by  charge  are 
implemented  as  arrays  of  floats  used  to  compute  the  aerodynamic  values. 

b.  Dynamic  parameters 

The  remaining  quantities  are  functions  of  mach  number,  and  tlie  intervals 
over  which  their  determining  coefficients  are  valid  must  be  obtained  differently  for  each 
one.  Therefore,  a  dynatnic  structure  is  used  to  contaui  the  sets  of  coefficients.  An  ordered 
list  of  objects  of  type  Physical_coefficient  is  created  for  each  parameter,  using  a  library  of 
container  classes  that  are  due  to  Wilson  (Wilson,  1 992).  The  list  is  ordered  by  the  maximum 
mach  number  for  which  that  particular  set  of  coefficients  is  valid. 

6.  Algorithm. 

a.  Initialization 

The  first  step  in  calculating  the  trajectories  is  initializing  the  dynamic  data 
structures.  Instances  of  objects  of  type  Physical_co€fficient  are  defined  that  correspond  to 
a  single  set  of  coefficients  needed  to  calculate  a  particular  aerodynamic  parameter  over  a 
certain  mach  interval.  The  values  for  these  coefficients  aie  taken  from  (Ballistic  Research 
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Laboratory,  1983,  pp.  2-8 ).  These  objects  are  inserted  into  a  list  in  order  of  the  highest  mach 
number  over  which  the  set  is  valid. 

b.  Numerical  integration 

The  remaining  step  of  the  algoritlnn  is  to  perform  numerical  integration  of 
Equation  3.8  over  the  interval  of  uiterest.  In  general,  this  will  be  from  r  =  0  until  impact. 
This  is  accomplished  using  a  predictor-corrector  technique  described  in  (Lieske  1973.  p. 
16).  A  predicted  value  iotf(i)  is  obtained  by  linear  extrapolation  from /(;-/)  by  applying  the 
first  derivative  over  the  width  of  the  integration  time  step.  This  estimate  is  used  to  calculate 
/'(/)  in  a  similar  fashion.  These  two  values  are  then  used  to  calculate  /  "(/),  which  is 
Equation  3.8.  The  result  is  used  to  calculate  corrected  values  of  the  first  two  derivatives  by 
the  addition  of  error-localizing  tenns  that  are  guaranteed  to  diminish  within  each  iteration 
of  the  computation.  The  values  thus  corrected  become  the  next  data  values  used  to  drive  the 
prediction  phase,  and  the  algorithm  continues  until  the  end  of  the  interval  is  reached.  S’mee 
the  ith  result  makes  use  of  values  from  the  /-/th  time  step,  it  is  necessary  to  retain  the  most 
recently  completed  values  from  each  iteration.The  instantaneous  values  for/t  are  obtained 
directly  by  using  Equation  3.15. 

The  values  for  the  aerodynamic  inputs  in  the  above  calculation  are 
determined  by  traversing  the  lists  that  were  created  as  described  in  Initialization  When  the 
set  is  found  that  is  valid  over  the  necessary  mach  interval,  the  value  for  the  parameter  is 
determined  by  the  following  equation; 

Pa  ram  =  aQ  +  a^M  +  a2M^  +  ...  +  (Eq3.16) 

where  a„  is  the  /ith  coefficient  of  the  set.  M  is  the  mach  number,  and  n  is  the  number  of 
coefficients  that  make  up  the  set. 
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A  comparison  of  the  results  obtained  using  this  variant  of  the  modified  point- 
mass  mod?’  is  shown  in  Table  2. 

TABLE  2;  MODIFIED  POINT-MASS  VARIANT  DATA  VS.  RRING  TABLE  DATA 


Charge 

Quadrant 

elevation 

(mils) 

Observed 

impact 

(deflection, 

range) 

Firing  table 
(deflection, 
range) 

CEP 

(meters) 

Error 

(number  of 
CEPs) 

3W 

266.67 

(26, 4077) 

(22,  4056) 

10.66 

2.0 

3W 

533.33 

(92.6555) 

(84.5,  6588) 

17.49 

1.93 

3W 

800.0 

(184. 7299) 

(207,  7300) 

21.35 

1.08 

7W 

266.67 

(63. 8940) 

(69. 8983) 

9.8 

4.43 

7W 

533.33 

(210. 

(220. 

14.97 

4.72 

13084) 

13154) 

7W 

800.0 

(453, 

(441, 

16.73 

2.67 

14657) 

14700) 
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IV.  TERMINAL  BALLISTICS 


A.  DEFINITION 

Terminal  ballistics  is  the  study  of  the  events  that  take  place  at  the  end  of  a  ballistic 
projectile’s  trajectory.  In  the  case  of  military  munitions,  the  desired  tenninal  effect  is 
usually  to  either  impart  an  enonnous  amount  of  kinetic  energy  to  the  target,  or  to  produce 
an  explosion. 

Because  of  the  difficulty  of  directing  them  with  precision  over  extended  distances, 
kinetic  energy  weapons  are  mainly  useful  against  low-signature,  heavily  armored  targets. 
The  total  energy  that  a  particular  projectile  could  theoretically  transfer  to  a  target  is  given 

1  ■> 

by  the  ubiquitous  formula,  e  =  -^m\^ .  Tire  means  through  which  these  weapons  function 

then  is  by  accelerating  a  very  dense  object  to  extremely  high  velocity. 

The  majority  of  ballistic  weapons  therefore  have  been  developed  for  use  against  area 
targets.  Upon  arrival  in  the  target  area,  these  munitions  are  designed  to  produce  lethal 
effects  by  means  of  the  rapid  combustion  of  volatile  compounds. 

B,  CHARACTERISTICS  OF  EXPLOSIVE  MUNITIONS 

Explosive  munitions  all  possess  a  common  attribute:  tliey  produce  destruction  in  the 
target  area  by  converting  then  stored  potential  energy  into  kinetic  energy  in  the  form  of 
blast,  heat,  and  light.  They  are  further  categorized  by  the  manner  in  which  this  conversion 
is  controlled. 

1.  Chemical  energy  weapons. 

Some  projectiles  are  designed  such  that  upon  impact,  an  intense  jet  stream  of  hot 
gases  is  projected  against  the  target.  Tlie  enonnously  high  temperatures  of  die  gases  cause 
portions  of  the  target  to  vaporize,  allowing  the  jet  to  penetrate  further,  which  causes  damage 
to  the  interior  of  the  target  area.  In  order  to  function  in  this  manner,  the  explosive  elements 
of  the  projectile  must  be  fonned  into  a  parabolic  shape,  which  focuses  the  energy  converted 
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by  the  explosion  into  the  characteristic  jet  stream  oriented  toward  the  target.  Thus,  these 
explosives  are  known  as  shaped  charge  munitions. 

Because  shaped  charge  munitions  depend  upon  concentrating  their  energy 
against  a  particular  place  on  the  target,  they  are  chiefly  useful  against  point  targets,  such  as 
bunkers  and  armored  vehicles. 

2.  Shrapnel-producing  weapons. 

The  majority  of  the  destructive  effect  in  tliese  kinds  of  weapons  is  attributable  to 
the  blast  component.  The  volatile  elements  are  surrounded  by  some  form  of  solid  matter 
that  disintegrates  into  smaller  pieces  under  the  forces  of  the  rapidly  expanding  gases 
produced  by  the  explosion.  The  kinetic  energy  of  the  gases'  expansion  is  imparted  to  the 
particles  of  sluapnel,  which  are  propelled  from  the  locus  of  the  explosion  at  high  velocity. 
The  shrapnel  then  collides  with  objects  in  the  target  area,  producing  the  desired  effect. 

From  the  virtual  world  perspective,  the  shrapnel-producing  weapons  pose  the 
more  interesting  problem:  how  can  the  visual  characteristics  of  such  weapons  be  modeled? 
To  answer  this  question,  we  must  examine  the  results  of  experimental  testing. 

Experiments  have  been  conducted  wherein  a  patticular  munition  is  detonated, 
and  the  final  positions  of  the  resultant  particles  of  sluapnel  are  recorded.  The  evidence 
shows  that  the  characteristics  of  these  particles  vary  with  the  shape  of  the  projectile,  the 
nature  of  the  confining  material,  and  the  type  of  explosive  (Reche,  1980,  pp.  170-171).  In 
general,  though,  the  initial  velocity  vectors  of  these  particles  can  be  assumed  to  have  a 
normal  distribution  with  respect  to  their  magnitudes,  and  a  uniform  distribution  with 
respect  to  their  orientations.  The  sizes  of  the  particles  produced  can  be  assumed  to  be 
normally  distributed.  Assuming  these  orderings  then,  a  key  element  in  the  development  of 
a  way  to  model  a  munition's  terminal  effects  in  a  synthetic  environment  is  appropriate 
selection  of  the  mean  and  variance  of  each  distribution. 

Another  important  attribute  in  the  visual  modeling  of  exploding  slirapnel  has 
been  referred  to  as  the  directionalit)’.  This  term  is  used  to  describe  the  global  orientation  of 


the  particles.  It  is  defined  by  a  principle  orientation  and  a  limiting  angular  deflection  known 
as  the  ejection  angle  (Reeves.  1983,  p.  367).  As  in  the  case  of  velocities,  a  normal 
distribution  of  shrapnel  pieces  within  the  ejection  angle  is  expected.  For  ground  bursts,  the 
surface  upon  which  the  munition  impacts  usually  forces  the  principal  orientation  to  be 
perpendicular  to  it.  Figure  2  shows  a  schematic  of  these  fundamental  quantities. 


Figure  2;  Essential  Visual  Cliaiacleristics  of  a  Shrapnel-Producing  Explosion  (Ground  Burst). 

In  the  case  of  air  bursts  the  directionality  is  neutral,  since  shrapnel  is  ejected 
roughly  equally  in  all  directions. 

C.  EXPLOSIONS  EDITOR  (NPSEE) 

As  an  aid  to  selecting  values  for  the  variables  described  in  the  previous  paragraphs, 
the  Naval  Postgraduate  School  Explosions  Editor  (NPSEE)  was  developed.  This  is  an 
interactive  design  tool  which  provides  a  means  to  map  out  the  overall  characteristics  of  an 
explosion  animation  sequence.  In  addition  to  the  global  attributes  already  discussed. 
NPSEE  allows  a  designer  to  specify; 

•  Shape  of  the  region  in  which  particles  aie  initially  distributed. 

•  Density  of  the  initial  particle  distribution. 
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•  Nature  of  the  particles  ’  individual  trajectories  (linear  or  parabolic ), 

•  Initial  color  of  the  particles. 

As  the  values  for  the  parameters  of  interest  are  adjusted,  an  on-screen  “template"  is 
changed  to  give  the  user  visual  feedback  of  the  effect  of  a  particular  action.  The  parameters 
are  then  saved  to  a  file,  which  can  be  read  in  to  control  the  actual  animation.  The  editor  is 
not  a  keyframing  or  sequencing  device.  It  does  not  seek  to  define  exactly  what  the  animated 
effect  will  look  like;  rather,  it  is  a  support  tool  that  can  be  used  to  test  various  combinations 
of  the  controlling  variables. 

1.  Interface. 

For  the  most  part,  input  to  NPSEE  is  managed  through  use  of  a  screen-oriented 
graphical  user  interface,  called  NFS  Pajiel  Designer  (NPSPD)  (King  &  Prevatt,  1990). 
NPSPD  produces  an  interface  that  is  intuitive  and  visual.  Users  provide  input  by 
manipulating  various  controls  on-screen  with  the  mouse.  The  controls  used  in  NPSEE  are 
either  buttons  or  sliders.  Pressing  the  left  mouse  button  while  the  cursor  is  over  a  slider  and 
then  dragging  the  slider  to  the  desired  position  allows  values  along  a  continuous  spectrum 
to  be  specified.  Some  controls  pertain  to  parameters  that  have  no  easily  implemented  visual 
analogue,  and  thus  changing  them  does  not  affect  the  shape  that  is  displayed  in  the  template 
window.  They  are  nonetheless  useful,  in  that  the  pomt-and-click  method  of  input 
management  is  usually  less  tedious  and  error-prone  than  typing.  Figure  3  shows  the  user 
input  screen  of  NPSEE. 

2.  Parameter  encoding. 

Visual  feedback  concerning  the  effect  of  changes  to  the  global  state  variables  is 
provided  by  means  of  an  iconic  shape  in  the  template  window.  It  is  intended  to  present  tire 
user  with  an  encapsulation  of  the  effects  that  the  current  parameter  settings  will  have  on  the 
explosion  particle  system.  This  is  accomplished  by  encoding  key  information  in  the 
characteristics  of  the  template  shapes  that  are  displayed.  For  an  air  burst,  the  particles  are 
flung  out  in  all  directions.  Thus,  a  spherical  shape  is  used  to  convey  this  type  of  dispersion. 
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In  a  ground  burst,  the  particles  are  mainly  propelled  outward  and  upward,  away  from  the 
site  of  the  explosion.  A  more  complex  shape  is  needed  to  capture  the  information  in  this 
case.  For  circular  ground  bursts,  a  roughly  cylindrical  shape  similar  to  that  shown  in  Figure 
3  is  constructed  using  non-uniform  rational  b-splines  (NURBS).  The  resulting  form 
conveys  information  in  the  following  manner: 

•  The  height  of  the  shape  determines  the  mean  magnitude  of  the  particles’  velocity 
distribution.  It  is  adjusted  using  the  ‘Y-mag’  (Y-magnitude)  button. 

•  The  angle  that  the  sides  of  the  shape  make  with  the  x-z  plane  determines  the  degree  to 
which  particles  are  cast  outward  from  the  center  of  the  explosion.  It  is  adjusted  using 
the  ‘Ejection  angle’  slider. 

•  The  overall  scale  of  the  shape  determines  the  size  of  an  individual  particle  for  both 
ground  and  air  burst.  It  is  adjusted  using  the  ‘XZ-mag’  slider. 


Figure  3:  The  NPSEE  Display. 

For  the  ground  burst  template,  three  viewing  options  are  available  by  pressing  the 
right  mouse  button  to  view  the  main  menu.  Either  the  splines  that  determine  the  shape,  the 
polygons  used  to  tessellate  the  surface,  or  the  interpolating  surface  itself  may  be  shown. 
Once  the  variables  have  been  adjusted  as  desired,  they  may  be  saved  to  a  file  for  use  in  the 
animation  sequence. 
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D.  PARTICLE  SYSTEMS 


1.  Fuzzy  objects. 

There  exist  physical  phenomena  occurring  in  nature  that  are  composed  of  a  large 
number  of  bodies  which  are  small  m  comparison  to  the  space  occupied  by  the  sum  of  the 
bodies  taken  as  a  whole.  Some  examples  are  clouds,  hair,  smoke,  fog,  steam,  fire,  dust,  and 
snow.  The  visible  manifestation  of  these  objects  have  a  differing  cross-sectional  density, 
usually  becoming  more  dense  toward  the  interior  of  the  volume  that  the  object  occupies. 
The  diminishing  of  body  density  toward  the  outer  limits  of  the  object  has  the  effect  of 
causing  the  edges  to  be  ill-defmed. 

Rendering  this  class  of  objects  in  a  virtual  world  poses  special  problems.  Tire 
computational  burden  of  determining  the  motion  of  each  body  during  the  screen  refresh 
period  is  too  great,  given  the  state  of  the  art  of  computer  hardware.  This  difficulty  is 
exacerbated  by  objects  whose  component  bodies’  visual  characteristics  also  change  over 
time.  In  general,  this  category  of  objects  incorporate  too  much  detail  to  represent 
convincingly  using  deterministic  methods. 

Aside  from  the  conunon  tltread  of  enormous  visual  detail,  most  of  the  entities  that 
comprise  these  “fuzzy”  objects  have  a  dynamic  lifetime,  the  period  of  time  throughout 
which  the  particles  are  visible.  As  an  example,  consider  a  cloud  of  smoke.  Over  time,  the 
cloud  moves  about,  and  is  dissipated  by  the  motion  of  the  air,  as  well  as  other  influences. 
When  a  portion  of  the  cloud  has  thinned  past  a  certain  density,  that  portion  appears  to  an 
observer  to  have  faded  away.  The  observer  will  notice  that  in  general,  clouds  of  smoke  tend 
not  to  fade  away  at  the  same  rate.  The  same  thing  is  true  of  fireworks.  Thus  there  is  a  need 
to  attenuate  a  particle’s  participation  in  the  simulation  over  time. 

2.  Management  of  complexity. 

A  particle  system  is  an  abstract  computational  structure  that  encapsulates  some 
of  the  complexity  of  this  class  of  objects.  They  are  collections  of  either  other  systems,  or 
atomic  particles.  They  necessarily  incorporate  some  form  of  .stochastic  algorithm  for 
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managing  the  aspects  of  the  object’s  appearance  and  behavior  that  cannot  be  accomplislieil 
explicitly.  The  attributes  that  are  frequently  governed  in  this  way  are  position,  velocity, 
size,  color,  shape,  and  Ufetime.  In  cases  where  the  object  is  an  aggregate  of  particle 
systems,  it  is  often  convenient  to  establish  a  hierarchical  relationship  between  them.  In  this 
way,  some  of  the  subordinate  elements’  global  qualities  can  be  encoded  into  their  positions 
in  the  hierarchy.  This  simplifies  the  task  of  managing  attributes  that  change  over  time. 

Exactly  how  to  determine  the  state  of  such  characteristics  as  a  function  of  time  is 
closely  tied  to  the  physical  phenomenon  being  modeled.  It  may  be  possible  to  discover  an 
analytic  function  that  will  suffice.  An  example  of  this  tecluiique  would  be  use  of  equations 
3.1  -  3.3  to  compute  the  position  of  a  projectile.  AJtematively.  a  linear  interpolation  might 
be  useful,  as  in 

/ig  =/(^„-l)  +  (Eq4.I7) 

Most  problem  domains  require  some  degree  of  empirical  evaluation  to  determine 
the  most  appropriate  technique.  In  situations  where  a  probabilistic  distribution  of  values  is 
required,  a  variant  of  the  technique  given  by  Reeves  (Reeves,  1983,  p.  361)  may  be 
suitable; 

f{t)  =  {RandQ  xa<j^)  (Eq4.l8) 

and  Of  are  the  mean  and  variance  of  die  desired  distribution,  respectively. 
Rand{)  is  a  function  that  returns  some  floating-point  value  between  0  and  1.  Clearly,  it  is 
important  that  the  implementation  of  Rand( )  be  chosen  to  reflect  any  modalities  that  might 
be  present  in  the  way  that  the  attribute  varies  in  reality.  The  method  of  constraining /(r)  to 

lie  within  ±ao  of  p  is  useful  when  fit)  models  a  physical  process  that  has  predictable 
boundary  conditions. 
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E.  IMPLEMENTATION 

The  notion  of  a  class  in  C++  is  an  abstraction  that  is  well-suited  to  the  implementation 
of  an  explosion  particle  system  model.  Tlie  ability  to  encapsulate  data  elements,  as  well  as 
defining  initialization  parameters  within  object  constructors  turns  out  to  have  direct 
analogues  to  the  particle  system  metaphor.  By  conceptually  equating  a  particle  with  an 
object,  the  logic  of  the  main  routine  that  accomplishes  the  motion  of  the  particles  is  greatly 
simplified. 

The  following  attributes  are  handled  in  a  stochastic  fashion: 

•  Initial  positions. 

•  Muzzle  velocities. 

•  Orientation  and  quadrant  elevation  used  to  calculate  the  particles’  trajectories. 

The  locations  of  the  particles  me  determined  such  that  they  lie  initially  within  the 
boundaries  of  a  circular  region.  The  formula  given  by  Equation  4.18  is  used  to  generate  a 
placement  set  that  is  square;  these  values  are  subsequently  clipped  to  lie  within  a  circle  of 
an  appropriate  radius.  The  distribution  of  these  locations  is  roughly  unifonn  with  respect  to 
the  circle’s  area. 

I.  Visual  characteristics. 

Each  individual  particle  in  the  simulation  is  a  tetrahedron.  The  NPSOFF  utilities 
discussed  above  were  used  to  render  these  objects.  Tlie  original  choice  was  a  two- 
dimensional  circle,  chosen  for  its  simplicity  and  consequent  speed  of  rendering.  The 
resulting  animation,  while  very  fast  in  terms  of  frames  per  second,  was  not  very  persuasive. 
To  make  the  picture  interesting,  it  seemed  that  a  form  that  had  shape  in  three  dimensions 
was  required.  The  next  logical  choice  was  a  polygonalized  sphere.  The  sphere  was  fairly 
expensive  to  render,  and  it  also  had  the  disadvantage  of  being  unbelievable  as  a  form  meant 
to  represent  a  piece  of  shrapnel.  The  tetralietlron  is  a  reasonable  compromise  between  these 
concerns. 
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2.  Data  structures. 


In  this  simulation,  a  system  of  particles  is  implemented  as  a  class  containing  a 
doubly-linked  list  of  atomic  objects.  Once  again,  good  use  is  made  of  Wilson's  NPSCL  to 
obtain  this  functionality,  and  to  maintain  the  object-oriented  motif  of  the  program.  Three 
objects  of  this  class  are  instantiated:  one  holds  the  state  of  the  system  as  it  is  being  rendered 
by  the  graphics  hardware;  one  holds  the  state  of  the  system  as  updates  are  calculated  for  the 
succeeding  frame;  the  last  is  used  as  a  buffer  between  the  rendering  and  calculation  states. 

3.  Algorithm. 

Upon  start-up,  each  particle  in  the  system  is  assigned  characteristics  of  position, 
velocity,  lifetime,  and  quadrant  elevation.  Using  a  simplified  dynamics  model,  the 
projected  location  of  each  particle  is  calculated  instantaneously  based  on  its  current 
location  and  velocity  vector,  and  rendered  accordingly.  The  temporal  resolution  of  this 
calculation  was  chosen  so  as  to  minimize  jerkiness,  but  still  give  a  faithful  unpression  of  a 
ballistic  trajectory.  On  an  Iris  4D/240  "VOX,  particle  counts  of  up  to  approximately  one 
hundred  result  in  satisfactorily  realistic  motion.  Every  particle  acts  under  a  force  designed 
to  simulate  the  earth’s  gravity,  as  do  particles  of  slirapnel  in  a  real  explosion.  Collision  with 
the  ground  is  detected  by  examining  the  y-coniponent  of  the  particle’s  position.  An  elastic 
point-to-point  collision  is  simulated,  with  a  linear  damping  of  the  original  particle  energy. 

Parallelism  is  exploited  by  partitioning  the  program  into  two  MIMD  processes. 
One  process  (the  “producer”)  is  tasked  to  calculate  succeeding  positions  for  each  particle. 
The  remaining  process  (the  consumer)  reads  the  output  from  the  producer,  and  displays 
particles  at  the  corresponding  positions.  A  two-buffer  scheme  due  to  Dr.  Michael  Zyda  is 
modified  to  incorporate  the  dependency  upon  previous  particle  velocities,  and  used  to 
enforce  mutual  exclusion  between  the  shared  data.  Some  preliminary  testing  showed  that 
the  program’s  performance  could  probably  be  improved  by  managing  a  queue  of  buffers, 
albeit  with  attendant  added  complexity.  Figures  4  and  5  show  the  resulting  output. 
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Figure  4:  Initial  Particle  System  Disposition  (t  =  0.0,  n  =  100). 


Figure  5:  Development  at  t  -  2.14. 
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V.  THE  AUTONOMOUS  FORWARD  OBSERVER 

A.  AUTONOMOUS  AGENTS 
1.  General. 

The  creation  of  a  computer  program  that  successfully  imitates  the  human 
capacity  for  perception  and  thought  has  been  the  holy  grail  that  computet  scientists  in  the 
field  of  artificial  intelligence  have  pursued  since  the  beginning  of  the  discipline.  According 
to  some  at  least,  little  real  progress  has  been  made  toward  that  end  (.Wilkes,  1992,  p.  17). 
Certainly  inroads  have  been  made  in  the  intensely  studied  area  of  computer  vision.  Edge 
detection,  region-growing,  and  stereoscopy  are  now  well-developed  methods  for  allowing 
a  hardware  and  software  system  to  “see".  Also,  work  in  the  fields  of  digital  signal 
processing  and  natural  language  understanding  has  enabled  computer  systems  to  “hear” 
and  respond  to  the  human  voice  (Luger  and  Stubblefield.  1989.  pp.  378-379).  It  would 
appear  that  the  harder  problem  has  been  defining  a  means  by  which  the  data  thus  received 
can  be  represented  and  reasoned  about.  Many  such  methods  have  been  proposed;  see 
(Jackson,  1990)  for  a  comprehensive  ueatment  of  the  subject  of  knowledge  representation. 
The  fact  that  no  single  method  has  been  shown  to  be  useful  in  the  general  case  of  reasoning 
over  many  wide  problem  domains  suggests  that  perhaps  the  definitive  work  in  this  area  has 
yet  to  be  undertaken. 

In  the  conduct  of  interactive  computer  sunulations  uivolving  the  combat  of 
ground  forces,  it  is  frequently  desirable  to  provide  a  capability  to  enable  entities  that  can 
automatically  emulate  the  actions  of  a  particular  group  of  participants  from  either  side.  This 
reduces  the  number  of  humans  required  to  conduct  the  simulation,  which  is  a  generally 
desirable  characteristic.  This  capability  is  especially  useful  vvhen  special  skills  are  required 
to  interact  in  the  roles  concerned.  To  this  end.  an  expert  system  could  be  constructed  with 
the  reasoning  abilities  necessary  to  simulate  the  cognitive  contributions  of  its  human 
counterpart. 
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Portraying  the  presence  of  autonomous  forces  in  a  virtual  world  is  related  to  the 
problem  of  creating  autonomous  robots  that  can  interact  with  human  beings,  but  with  an 
important  difference;  the  sensory  outputs  that  are  available  to  humans  as  well  as  the  inputs 
themselves  must  be  simulated.  How  to  simulate  the  various  cormnunications  channels,  as 
well  as  what  information  ought  to  be  provided  to  an  autonomous  entity  has  been  the  subject 
of  some  debate. 

2.  Global  V.S.  local  orientation. 

The  research  in  the  area  of  intelligence  simulation  in  autonomous  agents  has 
focused  mainly  on  the  development  of  efficient  algorithms  that  manipulate  an  internal 
symbolic  representation  of  the  world  knowledge  available  in  order  to  pursue  a  pre¬ 
determined  set  of  goals  (Maes.  1990,  p.  49).  Known  as  “traditional”  AI,  this  approach  has 
been  considered  by  some  to  be  problematic  for  certain  dynamic,  real-time  enviionments, 
especially  in  the  face  of  contradictory  inputs  (Pylyshyn,  1987).  An  alternative  paradigm  has 
been  developed  in  recent  years  which  emphasizes  the  tight  coupling  of  perception  to  action, 
without  the  need  for  deferring  a  response  until  the  envuoitment’s  global  state  can  be 
determined.  An  example  of  this  localized  orientation  is  the  subsumption  architecture 
(Brooks,  1989,  p.  692).  Research  with  autonomous  robots  has  shown  that  surprisingly 
complex  behaviors  can  appear  in  systems  that  are  programmed  to  react  to  sensory  inputs 
with  little  recourse  to  global  strategies.  There  is  no  reason  to  expect  that  the  same  would 
not  be  true  of  autonomous  agents  in  the  context  of  a  virtual  world. 

3.  Suitability  for  synthetic  environments. 

Autonomous  agents  designed  using  the  non-traditional  mode  have  several 
characteristics  that  are  desirable  in  terms  of  their  implementation  in  the  context  of  a  virtual 
world.  First,  the  presence  of  intelligent  agents  in  a  simulation  implies  an  increase  in 
program  overhead  having  computational  complexity  of  at  least  Oin).  due  to  the 
requirement  to  emulate  human  sensory  systems.  Each  agent  extant  in  such  a  system  must 
be  updated  periodically  regarding  what  etui  be  seen,  what  can  be  heard,  whether  or  not  it 
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has  collided  with  another  object,  etc.  It  is  therefore  necessary  to  minimize  communications 
to  and  from  agents,  if  graphical  representations  of  their  actions  are  to  be  rendered  in  real 
time.  Localized  architectures  allow  the  decision-making  modules  to  be  incorporated  at  a 
low  level,  obviating  the  need  for  passing  infonnation  received  from  the  sensory  modules  to 
a  global  planning  process.  Simulations  involving  the  participation  of  large  numbers  of 
agents  (human  or  otherwise)  can  rapidly  overwhelm  the  computing  capacity  of  even  the 
most  powerful  systems.  A  movement  toward  distributed  knowledge  and  function  seems 
inevitable  as  the  expectations  of  complex  behaviors  in  autonomous  agents  rises  with  the 
acceptance  and  employment  of  computer  simulations  in  general. 

The  object-oriented  paradigm  of  programming  is  widespread  in  the 
implementation  of  virtual  worlds.  Due  to  its  built-in  infonnation-hiduig  features,  the 
object-oriented  technique  is  a  natural  choice  to  implement  the  decentralized  decision¬ 
making  model.  Indeed,  it  is  assumed  that  complete  world  knowledge  is  not  available  to 
each  entity  in  such  an  enviromnent.  Not  only  does  this  have  the  effect  of  minimizing 
information  transmission,  it  encourages  the  development  of  a  standardized  interface  that 
allows  interaction  between  agents  that  occupy  a  simulation. 

As  graphical  simulations  become  more  sophisticated,  increasing  user 
expectations  of  realism  will  require  the  interaction  of  hundreds,  perhaps  thousands  of 
simultaneous  participants.  It  seems  unlikely  that  advances  in  computer  architecture  will 
result  in  the  creation  of  a  single  device  that  would  be  able  to  accommodate  the 
computational  load  that  such  a  simulation  would  entail,  at  least  not  in  real  time.  Therefore 
it  is  clear  that  any  simulation  that  seeks  to  maximize  the  appearance  of  realism  will 
necessarily  involve  the  cooperation  of  multiple  networked  computers.  This  suggests  that 
the  overall  nature  of  such  a  vSimulation  wiU  be  highly  parallel  and  decentralized.  The  non- 
traditional  approach  to  the  design  of  intelligent  agents  is  well  suited  to  such  an 
implementation. 


4.  Belief  systems. 

Another  design  decision  that  unplemeiitors  of  autonomous  entities  face  is  how  to 
determine  just  what  it  is  that  a  particular  agent  knows,  and  what  it  believes.  Note  that  in  this 
context  belief  is  the  agent’s  perception  of  fact  regardless  of  its  truth  in  the  absolute  sense; 
knowledge  is  awareness  of  absolute  truths  coupled  with  a  belief  to  that  effect.  To  some 
extent  the  manner  in  which  knowledge  and  beliefs  are  gained  is  related  to  the  choice  of 
orientation  discussed  above.  The  orientation  of  the  overall  decision-making  paradigm 
(global  planning  vs.  localized  reaction  to  stimulus)  will  affect  what  an  agent  is  permitted  to 
know  about  his  environment. 

When  a  goal  of  an  agent’s  participation  is  to  imitate  the  actions  of  some  human 
counterpart  (as  is  usually  the  case),  it  is  most  often  considered  undesirable  that  the  agent 
have  knowledge  about  the  world  that  is  complete  and  accurate.  This  is  tme,  of  course, 
because  human  beings  are  not  all-knowing.  Not  only  is  it  the  case  that  we  are  not  oracles, 
but  even  when  provided  with  explicit,  accurate  knowledge,  humans  cannot  always  be 
counted  upon  to  make  logical  decisions.  Tliis  is  attributable  to  a  variety  of  internal  and 
external  factors.  Most  of  the  external  factors  bearing  upon  Judgement  have  to  do  with  the 
attenuation  of  perception.  For  example,  a  soldier  that  is  wearing  night  vision  goggles 
(NVGs)  perceives  a  reconstruction  of  the  actual  world  unage  that  is  illuminated  in  such  a 
way  as  to  cause  the  appearance  that  every  object  is  some  shade  of  green.  Thus  his  sense  of 
vision  is  sharply  attenuated.  If  the  same  soldier,  still  wear  ing  NVGs,  was  driving  a  truck  in 
an  urban  area  where  there  were  color-coded  traffic  lights,  we  would  perhaps  expect  that  his 
imperfect  knowledge  about  the  color  of  the  lights  that  he  sees  might  lead  him  to  take  an 
inappropriate  course  of  action,  such  as  proceeding  through  an  intersection  when  the  traffic 
light  is  red.  (For  the  sake  of  illustration,  we  neglect  the  fact  that  an  additional  relational 
encoding  exists  with  respect  to  the  positions  of  the  lights  in  the  vertical  plane.) 

Thus  we  conclude  that  in  order  to  preserve  the  appearance  of  human  intelligence, 
the  sensory  apparatus  of  autonotnous  agents  must  somehow  be  distorted.  However,  not  just 
any  means  of  distortion  will  suffice.  It  is  a  safe  assumption  to  say  that  since  the  autonomous 


agents  we  are  discussing  are  meant  to  sunulate  the  actions  of  human  analogues,  then  in 
general  they  are  programmed  to  respond  to  sensory  stimuli  in  tlie  same  ways  that  human 
beings  do.  Therefore  it  is  not  reasonable,  for  example,  to  expect  a  program  that  controls  an 
agent’s  responses  to  consider  the  possibility  of  a  cow  that  it  encounters  producing  an 
automatic  weapon  and  taking  him  under  fire  (although  stranger  things  have  happened  in  the 
wide  world  of  virtual  realities  -  see  (Zyda  and  Pratt,  1992)).  In  iecting  a  filter,  it  must  be 
understood  that  too  large  a  deviation  from  tlie  actual  inputs  might  appear  obviously 
ludicrous,  and  destroy  the  illusion  of  realism.  Sensory  defonnations  that  give  the 
impression  that  a  law  of  physics  has  been  violated,  or  that  vastly  contradict  the  agent’s 
previous  experience  are  the  most  likely  to  produce  inappropriate  behavior  in  autonomous 
agents,  because  it  is  precisely  by  these  anticipated  events  that  its  responses  are  defined. 

Bhargava  and  Branley  have  explored  this  issue,  and  discovered  a  series  of  filters 
that  work  well  with  autonomous  agents  in  a  virtual  world  (Bhargava  and  Branley,  1992,  p. 
14).  Their  belief  function  is  designed  so  that  it  varies  logaritlunically,  and  may  be  clamped 
to  maximum  and  minimum  bounding  values  that  have  conespondence  to  the  appropriate 
limits  in  the  real  world.  They  also  describe  a  way  that  the  function  may  be  applied  so  as  to 
accommodate  contradictory  inputs,  usmg  Dempster-Shafer  theory  (Shafer,  1976).  Their 
method  was  used  to  distort  information  that  was  given  to  a  platoon  of  vehicles  acting 
autonomously  within  NPSNET,  with  very  satisfactory  results  (Branley,  1992). 

B.  AUTONOMY  IN  NPSNET 

All  vehicles  in  NPSNET  behave  autonomously,  with  the  exception  of  a  single 
“driven”  vehicle  per  workstation.  Tlie  diiveii  vehicle  can  be  changed  at  any  time  during  the 
simulation.  The  other  vehicles  are  programmed  to  meander  about  the  terrain  at  varying 
speeds  and  directions,  reversing  theix  diiection  when  either  the  virtual  ocean  or  the  end  of 
the  terrain  database  is  encountered.  Tliey  have  a  sensory  capability  that  makes  them  aware 
of  other  vehicles  that  are  firing  at  tliem.  iuid  can  interact  with  other  vehicles  by  either 
returning  fire  upon  them,  or  running  away  from  them. 
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Culpepper  implemented  a  module  in  NPSNET  that  allowed  for  a  platoon  of 
autonomous  armored  vehicles  to  operate  in  a  goal-oriented  fashion  (Culpepper.  1992). 
They  obtained  mission  and  subgoal  assignments  dynamically  in  response  to  changes  in  the 
environment,  such  as  the  presence  of  land  features,  and  incoming  weapons  fire.  These 
vehicles  also  have  the  capability  of  returning  fire,  as  well  as  initiating  offensive  firing  in  a 
direct  mode. 

The  autonomous  entities  discussed  so  far  are  limited  to  firing  upon  only  those  vehicles 
with  which  they  have  direct  visual  contact.  Tliis  is  by  conscious  design,  not  oversight;  the 
same  thing  is  tme  of  vehicles  operated  by  human  beings.  However,  there  is  a  way  that  fires 
can  be  brought  onto  objects  by  entities  beyond  the  targets'  line  of  sight. 

C.  DIRECT  VS.  INDIRECT  HRES 

As  in  real  combat,  participants  in  the  conduct  of  direct  fires  are  at  significant  risk  of 
being  acquired  as  targets  themselves  by  the  vehicle  that  is  beuig  fired  upon.  An  alternative 
to  enduring  this  risk  is  the  employment  of  indirect  fires,  as  the  field  artillery  does.  This 
technique  may  be  used  to  lob  munitions  onto  an  enemy  along  an  arcing  path,  rather  than 
propelling  them  straight  toward  the  target.  It  has  the  advantage  of  being  unhindered  by 
intervening  objects  such  as  terraui  features,  as  well  as  reducing  the  previously  mentioned 
risk  of  revealing  the  firing  element’s  location  to  the  enemy. 

Along  with  the  advantage  of  the  shooter’s  security  comes  an  equaUy  significant 
disadvantage:  because  of  the  characteristically  high  muzzle  velocity  of  military  weapons, 
the  range  to  target  in  general  is  beyond  the  capability  of  the  firing  entity’s  ability  to  observe. 
Thus,  it  is  usually  impossible  for  an  individual  using  a  weapon  that  fires  in  a  purely  indirect 
mode  to  acquire  targets  on  his  own  that  he  can  shoot  at  effectively.  As  well,  he  cannot 
evaluate  the  effects  of  the  fires  he  produces,  since  they  land  beyond  the  limits  of  his  vision. 
It  is  therefore  clear  that  in  order  to  employ  such  weapons  with  any  efficiency,  an  additional 
player  must  be  involved. 
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D.  THE  ARTILLERY  FORWARD  OBSERVER 

Tlie  person  responsible  for  acquiring  targets,  communicating  their  location  to  the 
firing  elements,  and  evaluating  the  effects  of  their  fires  against  the  targets  is  the  forward 
observer  (FO).  Since  the  FO  is  not  actually  firing,  it  is  possible  for  liim  to  remain 
undetected  while  attacking  the  targets.  In  order  to  perform  his  mission,  the  FO  must  be  able 
to  perform  the  following  cognitive  tasks; 

•  Be  aware  of  his  own  location. 

•  Sense  the  presence  of  other  entities  within  range  of  his  sensory  capabilities,  and 
determine  their  location. 

•  Determine  the  status  of  other  entities  as  either  friendly  or  enemy. 

•  Communicate  the  location  of  entities  determined  to  be  enemy  to  the  firing  elements. 

•  Sense  the  presence  of  munitions  arriving  iir  response  his  requests. 

•  Evaluate  the  effects  of  the  fires  noted,  and  take  appropriate  action  based  on  that 
evaluation. 

A  more  sophisticated  model  might  include  tasks  such  as  taking  actions  to  preserve  the 
observer’s  own  life,  should  the  need  arise.  Since  our  main  interest  is  in  reproducing  the 
capabilities  that  are  necessary  to  pennit  the  howitzers  to  fire  in  indirect  mode,  we  will 
neglect  such  issues. 

E.  IMPLEMENTATION 

All  of  the  tasks  listed  above  imply  the  receipt  or  transmission  of  sensory  infonnation 
as  a  prerequisite  to  the  reasoning  process.  Since  the  necessary  information  is  represented  in 
NPSNET  as  C  data  structures,  the  modules  that  simulate  the  sensory  processes  were 
implemented  in  the  C++  programming  language.  Once  the  sensory  infoimation  is 
processed  into  data,  it  is  communicated  to  an  appropriate  reasoning  module  written  in 
CLIPS.  These  modules  {rules,  in  CLIPS  parlance)  contain  the  logic  needed  to  perform  the 
observer’s  cognitive  tasks.  As  a  result  of  the  functioning  of  these  rules,  the  appearance  of 
an  autonomous  entity  that  behaves  in  the  manner  of  an  artillery  forward  observer  is 
obtained. 
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The  observer  is  made  aware  of  certain  key  default  data  at  program  start-up  by  means 
of  a  dejfacts  construct,  and  additionally  a  rule  that  creates  uistances  of  the  howitzers  and 
the  fire  direction  center.  The  deffacts  informs  the  observer  of  his  location,  and  the  kinds  of 
vehicles  he  is  to  consider  as  hostile.  The  howitzer  objects  are  made  aware  of  their  initial 
locations  as  well. 

The  fire  direction  center  (FDC)  is  an  entity  whose  pailicipation  in  real  indirect  tires  is 
crucial,  but  curiously  turned  out  to  be  of  minimal  significance  in  this  implementation.  In 
tactical  artillery  units,  persotuiel  in  the  fire  diiection  center  are  responsible  for  calculating 
the  orientation  of  the  howitzers  necessary  to  hit  the  target,  given  the  locations  of  the 
howitzers  and  the  target.  Since  this  is  a  very  mechanical,  algorithmic  process,  these 
calculations  are  performed  in  the  C-f-t-  modules,  and  the  FDC  object  is  maintained  in  the 
CLIPS  modules  mainly  for  the  sake  of  architectural  integrity  of  the  indirect  fire  paradigm. 

After  initialization,  NPSNET  has  a  main  driver  process  that  iteratively  detennincs 
changes  to  the  status  of  objects  (in  terms  of  their  location  and  appear  ance),  and  renders  the 
currently  visible  scene  on  the  user’s  display.  Code  is  inserted  into  this  process  such  that 
with  every  iteration: 

•  The  CLIPS  modules  are  advised  of  the  cunent  locations  of  the  howitzer  and  the 
observer. 

•  A  determination  is  made  of  what  vehicles  (if  any)  can  actually  be  seen  from  the 
observer's  position. 

•  The  CLIPS  modules  are  advised  of  the  locations  and  types  of  the  vehicles  that  can  be 
seen,  as  well  as  the  locations  of  impacting  artillery. 

•  The  CLIPS  inference  engine  is  allowed  to  run  to  completion. 

The  determination  referred  to  in  the  third  point  above  is  distinctly  non-trivial.  It 
involves  a  two-tiered  examination  of  every  vehicle  in  the  simulation.  First,  each  vehicle  is 
tested  to  see  if  it  is  close  enough  to  be  seen  with  the  naked  eye.  In  this  unplementation.  this 
distance  is  assumed  to  be  no  greater  than  four  kilometers.  All  vehicles  that  pass  this  first 
test  are  examined  for  the  possible  presence  of  obstructions  between  the  observer  and  the 
vehicle.  A  ray  is  traced  between  the  observer  and  the  target  by  determinuig  the  parametric 
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equations  of  a  line  between  those  two  points.  The  level  of  the  ground  is  calculated  at 
equidistant  intervals  along  each  ray  to  see  if  any  terrain  features  intervene.  The  locations 
and  types  of  the  vehicles  that  pass  both  tests  are  asserted  to  CLIPS  as  “contacts". 

The  CLIPS  modules  examine  each  contact  to  see  if  it  meets  the  criteria  to  be 
considered  a  target.  If  it  does  not,  the  contact  fact  is  simply  retracted.  If  the  contact  is  a 
hostile  one,  the  contact  fact  is  changed  to  a  target  fact.  Upon  assertion  of  the  target  fact, 
rules  determine  the  location  of  the  target  with  respect  to  tire  howitzers,  and  calculate  the 
necessary  deflections  and  elevations  needed  to  cause  the  rounds  to  land  as  close  to  the 
target  as  possible.  A  message  is  sent  to  the  howitzer  objects  to  slew  the  tubes  to  the  correct 
azimuth  and  elevation,  and  to  fire  themselves,  which  results  in  the  transmission  of  data  back 
to  NPSNET  to  do  so.  When  the  rounds  finally  land,  the  observer  sends  a  request  to 
NPSNET  to  determine  if  the  vehicle  that  it  was  shooting  at  was  in  fact  killed.  If  it  was 
killed,  then  the  observer  causes  the  guns  to  cease  firing.  If  it  was  not  killed,  the  observer 
causes  the  guns  to  fire  again.  If  the  target  is  in  motion  between  the  time  it  is  detected  and 
the  howitzer  fires,  the  observer  recalculates  the  mission  each  time  with  a  new  target 
location.  This  technique  differs  from  the  doctrine  currently  m  use  in  the  U.  S.  Army;  a  real 
observer  would  make  adjustments  in  tenns  of  deviations  from  the  locations  of  previous 
impacts,  rather  than  absolute  locations. 

F.  LIMITATIONS 

As  with  all  simulations  of  human  intelligence,  there  is  some  discrepancy  between  the 
operation  of  this  implementation  and  the  way  a  real  forward  observer  beliaves.  In 
particular; 

1.  Visibility  determination. 

The  ray-tracing  visibility  test  is  not  precise.  In  order  to  save  processing  time,  the 
elevation  of  the  ground  is  examined  at  no  more  than  one  hundred  locations  between  the 
observer  and  the  target.  Since  the  established  maximum  range  of  the  observer's  vision  is 
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approximately  four  thousand  meters,  this  means  that  terrain  features  of  less  than  forty 
meters  breadth  (measured  parallel  to  the  observer-target  line)  may  not  be  detected. 


2.  Range  resolution. 

The  increments  with  which  the  howitzers  can  adjust  their  elevation  is  limited  by 
the  length  of  the  lookup  table.  At  the  moment,  the  resolution  is  a  minimum  of  250  meters. 
This  makes  fast-moving  targets  difficult  to  kill.  A  numerical  root-finding  algoritlim  that 
solves  Equation  3.8  in  terms  of  quadrant  elevation  would  reduce  the  coarseness  of  this 
parameter. 

3.  Observer  behaviors. 

The  observer  is  limited  in  his  actions  upon  detection  of  a  target.  At  this  point,  the 
observer’s  operant  philosophy  could  be  summed  up  as.  “If  an  enemy  vehicle  is  detected, 
bring  fires  upon  it  until  it  is  dead.”  Real  observers  are  somewhat  more  flexible  in  their 
pursuit  of  the  mission.  While  the  format  of  the  facts  pertaining  to  mission  assignments 
requires  a  selection  from  one  of  the  three  classical  field  artillery  missions,  only  the  adjust 
fire  mission  is  supported  by  the  rule  base. 

The  observer  is  incapable  of  sensing  the  need  for  requesting  high-angle  fires,  as 
would  be  required  in  a  case  where  the  target  was  masked  from  the  firing  platoon  by  high 
terrain. 

4.  Battle  damage  assessment. 

The  kill  mechanism  of  a  bursting  projectile  is  represented  as  a  binary,  all-or- 
nothing  affair.  There  is  no  provision  for  assessing  levels  of  damage  less  than  complete 
destruction.  For  this  reason,  no  attempt  has  been  made  to  include  logic  that  selects  an 
optimum  shell-fuze  combination  based  on  the  target  chaiacteristics. 

5.  Belief  modes. 

The  sensory  data  received  by  the  observer,  and  consequently  by  the  FDC  to 
compute  the  firing  data,  is  perfect.  No  real-world  signals  attenuation  such  as  fog,  radio 


40 


static,  fatigue,  is  portrayed.  This  is  unrealistic,  because  many  of  the  subtasks  that  the 
observer  has  to  perform  depend  direcdy  on  the  quality  of  information  available  to  him. 
Examples  of  such  subtasks  are  target  acquisition  and  target  location.  Additionally,  there  is 
no  provision  for  possibilities  such  as  an  uiueliable  observer,  or  mistakes  in  judgement. 

Along  similar  lines,  an  assumption  here  is  that  no  misunderstandings  between  the 
FDC  and  the  howitzer  crews  take  place,  and  that  the  crew  orients  the  weapon  without  error. 
Experience  suggests  that  this  is  actually  an  appreciable  source  of  inaccuracy  in  weapons 
firing. 

These  limitations  notwithstanding,  the  implementation  lias  been  used 
successfully  in  a  small  combat  scenario  with  satisfactory  results.  Even  considering  the 
primitive  nature  of  the  autonomous  observer’s  abilities,  a  powerful  battlefield  force  can  be 
brought  to  bear  using  these  techniques  during  a  simulation. 
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VI.  CONCLUSION 


Finding  a  way  to  realistically  represent  objects  that  appear  in  a  particular  environment 
is  at  the  heart  of  all  simulation  systems.  Military  simulators  have  an  especially  acute  need 
to  find  representations  that  are  suitable  for  use  in  modeling  ballistic  projectiles,  since  these 
objects  are  always  present  in  non-trivial  combat  scenarios.  Simulations  that  seek  to  model 
military  operations  accurately  must  include  a  ballistics  model  that  is  as  faithful  to  the  actual 
weapon’s  performance  as  possible,  given  the  usual  speed  versus  accuracy  trade-offs. 
Lieske’s  modified  point-mass  trajectory  model  was  shown  to  be  suitable  for 
implementation  within  a  real-time  graphical  virtual  world  for  this  purpose. 

This  research  also  presents  a  means  by  which  the  visual  characteristics  of  the  terminal 
effects  of  shrapnel-producing  weapons  can  be  portrayed.  Although  the  particle  images 
themselves  lack  detail,  they  are  nonetheless  meaningful  because  the  richness  of  the 
variations  in  the  motion  of  the  particle  system  in  its  entirety  can  convey  an  animated 
sequence  that  is  believable  as  an  explosion.  This  can  enhance  the  realism  of  the  simulation 
overall.  Additional  research  is  necessary  to  discover  how  the  effects  of  the  destructive 
mechanisms  in  such  weapons  (i.e.,  energy  transfer)  can  be  similarly  modeled. 

The  field  artillery  is  known  in  some  circles  as  the  “King  of  Battle”,  allegedly  because 
it  is  responsible  for  the  highest  number  of  casualties  in  the  history  of  warfare.  In  any  event, 
the  ability  to  bring  fires  onto  a  target  area  that  is  miles  away  from  friendly  lines  is 
considered  fundamental  to  the  combat  tactics  of  every  modem  army  in  the  world.  The 
capacity  to  simulate  the  delivery  of  indirect  fires  is  therefore  an  important  part  of  combat 
simulators.  The  object-oriented  approach  toward  modeling  the  forward  observer  taken  in 
this  work  is  desirable  because  it  results  in  an  implementation  that  functions  in  much  the 
same  manner  as  its  human  counterpart.  This  produces  behavior  that  is  convincing,  through 
programming  that  is  easy  to  understand  and  modify. 
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APPENDIX  A  (VARIABLE  DEFINITIONS) 


A.  NOTATIONAL  CONVENTION. 

In  the  equations  used  in  this  paper,  vector  quantities  appear  in  bold  type,  with  a  vector 

symbol  directly  over  the  variable,  as  in  v.  All  other  variables  should  be  considered  to  be 
scalars.  When  a  variable  that  is  shown  in  non-boldface  type  without  the  vector  symbol  also 
appears  elsewhere  in  vector  format,  it  should  be  understood  that  the  instantaneous 
magnitude  of  the  vector  is  called  for.  This  avoids  the  somewhat  cumbersome  notation  of 

the  vector’s  norm,  as  in  ||  vl| . 

B.  VARIABLE  DEFINITIONS 

Equations  which  incorporate  variables  that  do  not  make  use  of  the  usual  SI  units  for 
that  quantity  (e.g.,  quadrant  elevation)  are  not  dependent  upon  them  for  their  dimensional 
correctness.  The  following  variables  are  used  throughout  this  paper: 


Variable 

Meaning 

Si.yni.ts, 

b 

Barrel  rifling  twist 

Calibers/rev 

C 

Ballistic  coefficient 

kg/m^ 

Drag  force  coefficient 

Dimensionless 

Co, 

Yaw  drag  force  coefficient 

Dimensionless 

Lift  force  coefficient 

Dimensionless 

% 

Yaw  lift  force  coefficient 

Dimensionless 

Spin  damping  moment  coefficient 

Dimensionless 

Cm. 

Overturning  moment  coefficient 

Dimensionless 

Cn,^ 

Magnus  force  coefficient 

Dimensionless 

d 

Diameter  of  the  projectile 

Meters 
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Meaning 


SI  units 


Variable 


g 

Gravitational  acceleration 
of  the  eanh 

Meters/ sec  ^ 

Axial  moment  of  inertia 

I  2 

kg  •  m 

L 

Lift  factor 

Dimensionless 

m 

Projectile  mass 

kg 

P 

Axial  spin  velocity 

Radians/second 

P 

Axial  spin  acceleration 

Radians /second 

Q 

Yaw  factor 

Dimensionless 

QE 

Quadrant  elevation 

Mils 

t 

Elapsed  time 

Seconds 

d 

Velocity  of  the  projectile  with 
respect  to  the  air 

Meters /second 

V 

Velocity  of  the  projectile  with 
respect  to  the  ground 

Meters/second 

«e 

Yaw  of  repose 

Radians 

A 

Coriolis  acceleration  of  the  earth 

meters /second^ 

P 

Air  density 

3 

kg  /  meter 
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APPENDIX  B  (USER’S  GUIDES) 


A.  BALLISTIC  TRAJECTORIES 

To  enable  the  calculation  of  ballistic  trajectories  in  NPSNET,  specify  the  “b”  option 
on  the  command  line,  for  example,  “npsnet  b”.  No  other  actions  are  required.  If  you  wish 
to  modify  the  default  muzzle  velocities  for  existing  weapons,  or  if  new  weapons  must  be 
added,  consult  the  file  ./data/weapons.dat.  It  is  an  ASCII  text  file  containing  the  parameters 
that  define  the  characteristics  of  weapons  in  NPSNET.  The  function  that  reads  this  file 
Ignores  comments  that  appear  which  use  the  “C”  programming  language  conventions.  The 
format  for  entries  in  the  file  is  given  here. 

The  identification  field  begins  with  the  character  “I”.  The  following  entries  must  be 
present,  in  the  following  order  and  format: 

•  char*  name  (name  of  the  weapon,  max  20  characters  with  no  spaces). 

•  int  id_num  (number  used  to  identify  a  particular  weapon  no  two  weapons  have  the 
same  id_num.  USSR  weapons  have  numbers  starting  from  2,  US  weapons  start  from 
100). 

The  data  field  begins  with  the  character  “D”.  The  following  entries  must  be  present, 
in  the  following  order  and  format: 

•  int  range  (absolute  maximum  range,  in  meters). 

•  float  muzzlevelocity  (how  fast  the  round  travels  as  it  leaves  the  firing  point,  in  meters 
per  second). 

•  int  killsize  (outside  radius  from  the  point  of  impact  of  lethal  effects,  in  meters). 

•  int  damagesize  (outside  radius  from  the  point  of  impact  of  non-  lethal,  but  nonetheless 
noticeable  effects  in  meters. 

•  char  knowledgesize  (outside  radius  from  point  of  impact  within  which  a  player  vehicle 
should  notice  that  a  round  has  impacted,  in  meters). 

•  char  targettype  (the  type  of  target  that  this  weapon  is  typically  most  effective  against. 
Key: 

~  “0”  =  Any  target 
~  “1”  =  Heavy  armor 
~  “2”  =  Light  armor 
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~  “3”  =  Wheeled  vehicle 
~  “4”  =  Personnel 
~  “5”  =  Aircraft 
~  “6”  =  Watercraft 

•  char  pathtype  (default  trajectory  employed  by  the  weapon.)  Key: 

~  “0”  =  low  angle  [direct  fire] 

-  “1”  =  high  angle  [indirect  fire] 

~  “2”  =  linear  [laser  beam] 

-  “3”  =  gravity  bomb. 

Note;  At  the  moment,  the  entries  for  killsize,  damagesize,  knowledgesize,  and 

targettype  are  not  used.  Several  examples  are  shown  below: 

I  1 5 SrratiMl 09 Howitzer  110  /*  Non  rocket-assisted.  Charge  3WB  */ 

D  7300  297.0  100  150  200  1  1 

I  2 03tranMl 10 Howitzer  111  /*8-inch  howitzer,  non  rocket  assisted*/ 
D  16800  370.0  150  200  250  1  1 

I  50cal  112 

D  1800  100.0  2  5  50  2  0 

IMPORTANT  NOTE:  The  autonomous  forward  observer  is  designed  to  operate  with 
the  firing  platoon  shooting  at  a  muzzle  velocity  corresponding  to  charge  three  white  bag. 
The  reason  is  that  the  algorithm  for  calculating  the  quadrant  elevation  for  the  howitzers  to 
achieve  the  range  necessary  to  hit  the  target  must  make  an  assumption  about  the  muzzle 
velocity  in  order  to  function,  given  that  the  achieved  range  is  a  direct  function  of  charge  and 
elevation.  If  the  default  muzzle  velocity  for  M 109  howitzers  is  changed  to  a  value  different 
than  charge  three,  the  howitzers  will  not  respond  properly. 

B.  NAVAL  POSTGRADUATE  SCHOOL  EXPLOSION  EDITOR 

The  meanings  of  the  graphical  controls  in  NPSEE  were  described  in  Chapter  IV  .  In 
addition,  the  following  keyboard  commands  are  in  use: 

•  Keyboard  up  arrow,  down  arrow;  viewpoint  z-,  z+. 

•  Keyboard  left  arrow,  right  arrow:  viewpoint  x+  x-. 
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•  Keyboard  home,  end:  viewpoint  y+,  y- 

•  Keypad  up  airow,  down  arrow  (8  and  2):  reference  point  y+,  y-. 

•  Keypad  left  arrow,  right  arrow  (4  and  6):  reference  point  x-,  x+. 

C.  AUTONOMOUS  FORWARD  OBSERVER 

To  enable  the  autonomous  forward  observer  in  NFS  NET,  specify  me  “a”  and  the  “b” 
options  on  the  command  line,  for  example,  “npsnet  a  b  t”.  No  other  actions  are  required. 
Upon  program  start-up,  the  forward  observer’s  vehicle  (an  M577  by  default),  and  a  platoon 
of  four  M109  howitzers  will  be  added  to  the  array  of  vehicles  that  is  either  read  in  from  a 
script,  or  initialized  randomly.  The  file  forward_observer.clp  must  be  modified  if  you  wish 
to  change  either  of  the  following  parameters: 

•  The  observer’s  initial  position,  and  the  vehicle  upon  which  he  is  mounted. 

•  The  initial  positions  of  the  howitzers  in  the  firing  platoon. 

•  The  types  of  vehicles  that  the  observer  is  to  recognize  as  hostile. 

The  location  of  the  observer  is  in  local  X-Z  coordinates.  This  vehicle  will  appear  in 
the  world  along  with  the  other  vehicles  that  are  read  in  from  ./data/vehposfiles/ 
vehiclepos.cnvy,  or  the  default  random  vehicles.  The  default  observer  location  is  the  center 
of  the  world,  as  in: 

(location  observer  25000.0  25000.0) 

The  kind  of  vehicle  that  the  observer  is  mounted  upon  is  an  index  into  the 
vehtypearray  data  structure  declared  in  ./headers/npsnet.h.  The  association  between 
indices  and  vehicles  can  be  found  in  the  file  ./datafiles/vehicle. dat.  The  default  is  ‘8’  for  an 
M577,  as  in: 

(mounted_on  8) 

The  observer’s  mission  is  specified  in  terms  of  vehicles  to  watch  for,  and  what  to  do 
when  he  sees  these  vehicles.  This  means  that  you  can  program  the  observer  for  stupidity. 
The  general  template  for  these  specifications  is: 

(mission  <action>  <vehicle-index>) 
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“<vehicle-index>“is  an  index  to  the  vehtypearray,  as  mentioned  above.  “<action>”  is 
either  ‘destroy’,  ‘neutralize’,  or  ‘special’. 

The  default  mission  is  just  to  kill  T-72s,  BMPs,  and  T2freds,  as  follows: 

(mission  destroy  51) 

(mission  destroy  52) 

(mission  destroy  53) 

The  location  of  each  gun  in  a  four-gun  firing  platoon  defaults  to  the  southeastern 
ponion  of  the  Fort  Hunter-Ligget  terrain  database: 

(location  gunl  16000.0  16000.0) 

(location  gun2  25505.0  25506.0) 

(location  gun3  24000.0  24000.0) 

(location  gun4  24500.0  24500.0) 
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